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Abstract 

We consider an elliptic self-adjoint first order differential operator act- 
ing on pairs (2-columns) of complex-valued half-densities over a connected 
compact 3-dimensional manifold without boundary. The principal sym- 
bol of our operator is assumed to be trace-free. We study the spectral 
function which is the sum of squares of Euclidean norms of eigenfunctions 
evaluated at a given point of the manifold, with summation carried out 
over all eigenvalues between zero and a positive A. We derive an explicit 
two-term asymptotic formula for the spectral function as A —5- +00, ex- 
pressing the second asymptotic coefficient via the trace of the subprincipal 
symbol and the geometric objects encoded within the principal symbol — 
metric, torsion of the teleparallel connection and topological charge. We 
then address the question: is our operator a massless Dirac operator on 
half-densities? We prove that it is a massless Dirac operator on half- 
densities if and only if the following two conditions are satisfied at every 
point of the manifold: a) the subprincipal symbol is proportional to the 
identity matrix and b) the second asymptotic coefficient of the spectral 
function is zero. 
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1 Main results 

Consider a first order differential operator A acting on 2-columns v = (y\ V2) T 
of complex- valued half-densities over a connected compact 3-dimensional man- 
ifold M without boundary. (See subsection 1.1.5 in for definition of half- 
density.) We assume the coefficients of the operator A to be infinitely smooth. 
We also assume that the operator A is formally self-adjoint (symmetric): 



w*Avdx= / (Aw)*vdx (1.1) 

M JM 

for all infinitely smooth v, w : M — > C 2 . Here and further on the superscript * 
in matrices, rows and columns indicates Hermitian conjugation in C 2 and dx := 
dx x dx 2 dx 3 , where x — (x 1 , x 2 , a; 3 ) are local coordinates on M. 

Let Ai(x, £) be the principal symbol of the operator A, i.e. matrix obtained 
by leaving in A only the leading (first order) derivatives and replacing each 
d/dx a by i£ a , a — 1, 2, 3. Here £ = (£i, £2, £3) is the variable dual to the position 
variable x; in physics literature the £ would be referred to as momentum. Our 
principal symbol A\(x 1 £) is a 2 x 2 Hermitian matrix-function on the cotangent 
bundle T*M, linear in every fibre T*M (i.e. linear in £). 

Throughout this paper we assume that the principal symbol A\ (x,^) is trace- 
free for all (x,0 e T*M and that 

detAi(x,0 t^O, V(i,0eT'M, (1.2) 

where T'M := T*M\{£ = 0} (cotangent bundle with the zero section removed). 
The assumption (|1.2[) is a version of the ellipticity condition. 

Under the above assumptions A is a self-adjoint operator in L 2 (M;C 2 ) 
(Hilbert space of square integrable complex-valued column "functions") with 
domain H l {M] C 2 ) (Sobolev space of complex- valued column "functions" which 
are square integrable together with their first partial derivatives) and the spec- 
trum of A is discrete, with eigenvalues accumulating to ±00. Let and 



2 



Vk = (vki(x) Vk2(x)) T be the eigenvalues and eigenfunctions of the opera- 
tor A. The eigenvalues A& are enumerated in increasing order with account of 
multiplicity, using a positive index k = 1, 2, . . . for positive A& and a nonpositive 
index k — 0, —1, —2, ... for nonpositive A&. 

We will be studying the spectral function and the counting junction. The 
spectral function is the real density defined as 

e(X,x,x):= IM*)!! 2 ' (1-3) 

0<A fc <A 

where ||vf.(a;)|| 2 :— [vk(x)]*Vk(x) is the square of the Euclidean norm of the 
eigenfunction Vk evaluated at the point x £ M and A is a positive parameter 
(spectral parameter). The counting function is the function 

N(X) := 1 = / e(X,x,x)dx. (1.4) 

0<A fc <A •' M 

In other words, N{\) is the number of eigenvalues Xk between zero and A. 

We aim to derive, under appropriate assumptions on Hamiltonian trajec- 
tories, two-term asymptotics for the spectral function (11.3[) and the counting 
function (|1.4[) . i.e. formulae of the type 

e(A, x, x) = a{x) A 3 + b(x) A 2 + o(A 2 ), (1.5) 

N(X) = aX 3 + bX 2 + o(X 2 ) (1.6) 

as A — > +cxd, where the real constants a, b and real densities a{x), b{x) are 
related in accordance with 

a(x)dx, (1-7) 

M 



b= \ b(x)dx. (1.8) 

In our recent paper |10| we performed a comprehensive analysis of two-term 
spectral asymptotics for general first order elliptic systems. In doing this we 
showed that all previous publications on systems gave formulae for the second 
asymptotic coefficient that were either incorrect or incomplete (i.e. an algorithm 
for the calculation of the second asymptotic coefficient rather than an explicit 
formula), see Section 11 of [10] for the appropriate bibliographic review. The 
correct formula for the coefficient b{x) was the main result of |10| . 

The problem examined in the current paper is a special case of that from [10J. 
Namely, in the current paper we make the following additional assumptions as 
compared to |10| : 

our manifold has dimension 3, (1-9) 
the number of equations in our system is 2, (1-10) 
our operator is differential (as opposed to pseudodifferential), (1-11) 
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the principal symbol is trace-free. (1-12) 

The need for a detailed analysis of the special case (|1.9|) - (|1.12p is driven by 
applications to the massless Dirac operator. 

The additional assumptions (|1.9I) - (|1.12[) lead to the following simplifications 
as compared to |10| . 

• The subprincipal symbol A sn b does not depend on the dual variable £ 
(momentum) and is a function of x (position) only. Recall that the sub- 
principal symbol is the zeroth order term of the full symbol of the first 
order operator A written in a way which makes it invariant under coordi- 
nate transformations, see formula (|6.2I) for formal definition and subsection 
2.1.3 in [55] for background material. 

• The principal symbol A\ admits a geometric description. 

The first of these simplifications is trivial whereas the second is not. We list 
below the geometric objects encoded within the principal symbol. 

Geometric object 1: the metric. Observe that the determinant of the 
principal symbol is a negative definite quadratic form 

detA 1 (x,0 = -g aP ^ (1-13) 

and the coefficients g a @(x) — g@ a (x), a,f3= 1,2,3, appearing in (| 1 . 1 3|) can 
be interpreted as components of a (contravariant) Riemannian metric. This 
implies, in particular, that our Hamiltonian (positive eigenvalue of the principal 
symbol) takes the form 

and the x-components of our Hamiltonian trajectories become geodesies. 

Geometric object 2: the teleparallel connection. This is an affine 
connection defined as follows. Suppose we have a covector £ based at the point 
x £ M and we want to construct a parallel covector £ based at the point x £ M. 
This is done by solving the linear system of equations 

A 1 (x,£)=A 1 (x,S). (1.15) 

Equation (| 1 . 1 5|) is equivalent to a system of three real linear algebraic equations 
for the three real unknowns, components of the covector £, and it is easy to 
see that this system has a unique solution. It is also easy to see that the 
affine connection defined by formula (|1.15[) preserves the Riemannian norm of 
covectors, i.e. g a ^{x)^ a ^p = g a @(x) £a£,3, hence, it is metric compatible. The 
parallel transport defined by formula (j!.15[) does not depend on the curve along 
which we transport the (co)vector, so our connection has zero curvature. The 
word "teleparallel" (parallel at a distance) is used in theoretical physics [18] to 
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describe metric compatible affiiie connections with zero curvature. The origins of 
this terminology go back to the works of A. Einstein and E. Cartan |28[ [23l [7]. 
though Cartan preferred to use the term "absolute parallelism" rather than 
"teleparallelism". 

The teleparallel connection coefficients T a p 7 (x) can be written down explic- 
itly in terms of the principal symbol, see formula (|3.7[) . and this allows us to 
define yet another geometric object — the torsion tensor 

T a j3 1 '-—T a ^ — Y a 1 p . (1-16) 

Further on we raise and lower indices of the torsion tensor using the metric. 

Geometric object 3: the topological charge. It turns out, see Sec- 
tion [31 that the existence of a principal symbol implies that our manifold M is 
parallelizable. Parallelizability implies orientability. Having chosen a particular 
orientation, we allow only changes of local coordinates x a , a = 1,2,3, which 
preserve orientation. 

We define the topological charge as 



c := --y/detg af ) tr((Ai) fl (A^ (A 1 ) Ss ) , (1.17) 

with the subscripts £ Q indicating partial derivatives. We show in Section [3] 
that the number c defined by formula (|1.17[) can take only two values, +1 or 
— 1, and describes the orientation of the principal symbol relative to the chosen 
orientation of local coordinates. 

We have identified three geometric objects encoded within the principal sym- 
bol — metric, teleparallel connection and topological charge. Consequently, one 
would expect the coefficient b(x) from formula (|1.5|) to be expressed via these 
three geometric objects and the subprincipal symbol. This assertion is confirmed 
by the following theorem. 

Theorem 1.1 The coefficients in the two-term asymptotics jl.5\) are given by 
the formulae 




(1.18) 



b(x) = — ([ 3 c *T ax - 2 tr A suh } ^det g afs ) (x) , (1.19) 

where 

T a0~f : = 3(^7 + T fa/3 + Tp la ) (1-20) 

is axial torsion (totally antisymmetric piece of the torsion tensor) and * is the 
Hodge star (3^ ). 

Remark 1.1 The spectral and counting functions admit two-term asymptotic 
expansions il.5\) and il.6]) only under appropriate assumptions on geodesic loops 



and closed geodesies respectively, see Theorems 8.3 and 8.4 in 'J Of. However, 
one can easily reformulate asymptotic formulae \1.5\) and (1.6\) in such a way 
that they remain valid without assumptions on geodesies: this can easily be 
achieved, say, by taking a convolution with a function from Schwartz space <S(R), 
see Theorems 7.1 and 7.2 in JTUj. Thus, the second asymptotic coefficients of 
the spectral and counting functions are well-defined irrespective of how many 
geodesic loops or closed geodesies we have. We introduced the second asymptotic 
coefficients b{x) and b via the unmollified asymptotic expansions il.5\) and il.6}) 
simply for the sake of clarity of presentation. 

The proof of Theorem 11.11 is given in Sections [5115] 

We now turn our attention to the massless Dirac operator. This operator is 
defined in Appendix see formula (IA.3[) . and it does not fit into our scheme 
because it is an operator acting on a 2-component complex- valued spinor (Weyl 
spinor) rather than a pair of complex- valued half-densities. However, on a par- 
allelizable manifold components of a spinor can be identified with half-densities. 
We call the resulting operator the massless Dirac operator on half- densities. 
The explicit formula for the massless Dirac operator on half-densities is (| A. 19|) . 

The massless Dirac operator on half-densities is an operator of the type we 
are considering in this paper, i.e. a self-adjoint first order elliptic differential 
operator acting on 2-columns of complex-valued half-densities and with a trace- 
free principal symbol. We address the question: is a given operator A a massless 
Dirac operator? The answer is given by the following theorem which is our main 
result. 

Theorem 1.2 The operator A is a massless Dirac operator on half-densities 
if and only if the following two conditions are satisfied at every point of the 
manifold M : a) the subprincipal symbol of the operator, A^i^x), is proportional 
to the identity matrix and b) the second asymptotic coefficient of the spectral 
function, b(x), is zero. 

Note that conditions a) and b) in Theorem 11.21 are invariant under special 
unitary transformations, i.e. transformations of the operator 

A^RAR*, (1.21) 

where R : M —> SU(2) is an arbitrary smooth special unitary matrix-function. 
The invariance of condition b) is obvious. In fact, condition b) is invariant under 
the action of a broader group: the unitary matrix-function R(x) appearing in 
formula (|1.2ip does not have to be special. As to condition a), its invariance 
is established by examination of formula (9.3) from [TU] with the use of the 
special commutation properties of trace-free Hermitian 2x2 matrices. The fact 
that the conditions of Theorem 11.21 are SU(2) invariant is not surprising as the 
massless Dirac operator is designed around the concept of SU(2) invariance, see 
Property 4 in Appendix [A] 

The proof of Theorem 11.21 is given in Sections [6] and [71 
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Theorems 11.11 and 11.21 tell us that for the massless Dirac operator on half- 
densities formulae (|1.5[) and (|1.6[) read 



e (^*) = ^A» + o(n (1.22) 

^( A ) = ^?- A3 + ( A2 )' ( L23 ) 

where VolM is the volume of the Riemannian 3- manifold M. 

Remark 1.2 The factor ydet g a p{x) appears in the RHS of U.22]) because we 
are working with the massless Dirac operator on half- densities HA. 19]) rather 
than with the massless Dirac operator on spinors HA.3\) . For the massless Dirac 
operator on spinors the spectral function is a scalar field (as opposed to a density) 
and formula U.22)) reads e(X,x,x) = g^-A 3 +o(A 2 ). 



2 Reduction from the general setting 

As explained in Section [1] the problem considered in the current paper is a 
special case of that from [TO]. Formulae (1.23) and (1-24) from [T0J in our case 
read 

a(x)= J d£, (2.1) 

fc+(x,f)<l 

b(x) = b 1 (x) + b 2 (x), (2.2) 

where 



6i(s)=-3 J ([v+}*A suh v+)(x,S)aX, (2.3) 
h+(x,£)<l 

k0«0 = ? [ {[v + ]*,A 1 -2h+I,v+}(x,0^. (2.4) 



h+(x,$)<l 

Here h + (x, £) is the positive eigenvalue of the principal symbol (see also formula 
p,14p ). v + (x,£,) is the corresponding normalized eigenvector (2-column), d£ is 
shorthand for d£ := (2tt)~ 3 d£ = {2tt)^ z d^id^d^z and / is the 2x2 identity 
matrix. Curly brackets in formula (jTTTSJ) denote the Poisson bracket on matrix- 
functions 

{P,R}:=P xa R Sa -P ia R xa (2.5) 
and its further generalization 

{P,Q,R}:=P xa QR u -P u QR xa , (2.6) 

with the subscripts x a and £ a indicating partial derivatives and the repeated 
tensor index a indicating summation over a = 1,2, 3. 
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Put P + (x,£) := [v + (x,l;)][v + (x,£,)]* , which is the orthogonal projection 
onto the eigenspace spanw + of the principal symbol. We have A\ — 2h + I = 
2h+P+ - Zh+I and {[«+]*, P+, v+} = 0, so formula can be rewritten as 

H*) = ~ I {h+{[v+]*,v+}){x,Z)dZ. (2.7) 



2 

h+(x,()<l 

Our aim now is to evaluate the integrals (|2.ip . (|2.3[) and (|2.7[) explicitly. 

Formulae (|2.ip and (|1.14p immediately imply (|1.18p . 

In order to evaluate the integral (|2.3p we rewrite this formula as 

h(x) = -3 J tv(A suh P+)(x,0^ 

and use the fact that P + (x,£) = 2h +^ x ^ (Ai(x,£) + h + (x,£) I). We get 
h(x) = -3 J 2/l+ L ^ tr(Amb(^i + h + I))(x, <g . 

h+(a: i e)<l 

But ^4 S ub does not depend on £ whereas Ax is an odd function of £, so the term 
tr(A su bAi) integrates to zero, leaving us with 

3 f 1 

h{x) = --(trA auh ){x) / <S£ = - — (tr A sub ^det ) (a?) . (2.8) 

h+(x,0<l 

In order to complete the proof of Theorem 1 1.1 1 we need to evaluate explicitly 
the integral (|2.7p . The next three sections deal with this nontrivial issue. 



3 Teleparallel connection 

We show in this section that the principal symbol generates a teleparallel con- 
nection which allows us to reformulate the results of our spectral analysis in a 
much clearer geometric language. 

Let us show first that the existence of a principal symbol implies that our 
manifold M is parallelizable. The principal symbol A\(x, £) is linear in £ so it 
can be written as 

A 1 (x,£) = a Q (x)£ Q , (3.1) 

where a a (x), a — 1,2,3, are some trace-free Hermitian 2x2 matrix-functions. 
Let us denote the elements of the matrices a a as a a ab , where the dotted index, 
running through the values 1,2, enumerates the rows and the undotted index, 
running through the values 1, 2, enumerates the columns; this notation is taken 
from [TT]. Put 

e 1 a {x) :=Rea a i2 (x), e 2 a {x) := -lma a i2 (x), e 3 a {x) := a a ^{x) . (3.2) 
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Formula (|3.2p defines a triple of smooth real vector fields ej(x), j = 1,2,3, on 
the manifold M. These vector fields are linearly independent at every point 
x of the manifold: this follows from formula (|1.2[) . Thus, the triple of vector 
fields ej is a frame. The existence of a frame means that the manifold M is 
parallelizable. 

Conversely, given a frame ej we uniquely recover the principal symbol Ai (a;, £) 
via formulae (|3.ip . (|A.1[) and (|A.2I) . Thus, a principal symbol is equivalent to a 
frame. Of course, this equivalence statement relies on our a priori assumptions 
(fTTH . (OP and (01)- (PTT2l . 

It is easy to see that the frame elements ej are orthonormal with respect to 
the metric (|1.13p . Moreover, the metric can be determined directly from the 
frame as 

g ^ = V k e« e k p , (3.3) 

where the repeated frame indices j and k indicate summation over j, k — 1, 2, 3. 
The two definitions of the metric, (|1.13l) and (|3.3[) . are equivalent. 

Parallelizability implies orientability, see Proposition 13.5 in [20J. Having 
chosen a particular orientation, we allow only changes of local coordinates x a , 
a = 1,2,3, which preserve orientation and define the Hodge star in the standard 
way: the action of * on a rank q antisymmetric tensor Q is 

(*Q) 7 , +1 ... 73 := W-r 1 v/^t^Q 71 - 7 ^..^ , (3-4) 

where e is the totally antisymmetric quantity, £123 := +1, and g is the Rie- 
mannian metric (|1.13p . Here and further on we identify differential forms with 
covariant antisymmetric tensors. We raise and lower tensor indices using our 
metric. 

Substituting formulae (|3.1|l and ()3.2p into (|1.17|) we get 

c = sgndet ej a . (3.5) 

Formula (|3.5p provides an equivalent (and more natural) definition of topological 
charge. It also explains why the topological charge, initially defined in SectionQ] 
in accordance with formula (|1.17p . can only take values +1 or — 1. 

The concept of a teleparallel connection was already defined in Section Q] in 
accordance with formula (|1 . 15|) . This connection can be equivalently defined 
via the frame as follows. Suppose we have a vector v based at the point x £ M 
and we want to construct a parallel vector v based at the point x € M. We 
decompose the vector v with respect to the frame at the point x, v = c J ej(x), and 
reassemble it with the same coefficients cP at the point x, defining v := cPej(x). 

We now define the covariant derivative corresponding to the teleparallel con- 
nection. Our teleparallel connection is a special case of an affine connection, so 
we are looking at a covariant derivative acting on vector /covector fields in the 
usual manner 

v, t v a = dv a /dx» + rv v p , = dwp/dx" - rv Wa . 
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The teleparallel connection coefficients are defined from the conditions 



V M e/*=0, (3.6) 

where the ej are elements of our frame. Formula (|3.6p gives a system of 27 linear 
algebraic equations for the determination of 27 unknown connection coefficients. 
It is known (see, for example, formula (A2) in [5]), that the unique solution of 
this system is 

T a tf = e k a (de k p/dx»), (3.7) 

where 

e k p := S ki 9 M ef. (3.8) 

The triple of covector fields e k , k = 1,2, 3, is called the coframe. The frame and 
coframe uniquely determine each other via the relation 

e 3 a e k a =8 k . (3.9) 

Note that our notation for the frame and coframe is taken from [T3] . We feel 
it necessary to mention this because there is whole range of different notation 
for frames/coframes in mathematics and theoretical physics literature, which 
makes the subject somewhat confusing. 

One can check by performing explicit calculations that the teleparallel con- 
nection has the following two important properties: 

V a g Pl = 0, (3.10) 

which means that the connection is metric compatible, and 

(Vq.V/3 — V^Vq,)^ 7 = for any vector field v , (3-11) 

which means that the Riemann curvature tensor is zero. Properties (|3.10j) and 
(|3.1ip are the defining properties of a teleparallel connection: a teleparallel 
connection is, by definition [18], an affine connection satisfying (|3.10p and (|3. 1 1|) . 

The tensor characterizing the "strength" of the teleparallel connection is not 
the Riemann curvature tensor but the torsion tensor (|1.16[) . The teleparallel 
connection is, in a sense, the opposite of the more common Levi-Civita con- 
nection: the Levi-Civita connection has zero torsion but nonzero curvature, 
whereas the teleparallel connection has nonzero torsion but zero curvature. In 
our paper we distinguish these two affine connections by using different notation 
for connection coefficients: we write the teleparallel connection coefficients as 
T a p-y and the Levi-Civita connection coefficients (Christoffel symbols) as | ^ \, 
see formula (|A.4[) . It is known, see formula (7.34) in [21], that the two sets of 
connection coefficients are related as T a p^ = | ^ j + i(T a ( g 7 + Tp a 1 + T 7 Q | g). 

Substituting (|3.7[) into (|1.16p we arrive at the following explicit formula for 
the torsion tensor of the teleparallel connection 

T^ej^de 1 , (3.12) 
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where the d stands for the exterior derivative. For the sake of clarity we rewrite 
formula (|3 . 1 2|) in more detailed form, retaining all tensor indices, 

T a M = e^ide^/dx? - de j /dx^) . (3.13) 

As always, the repeated index j appearing in formulae (|3.12|) and (|3.13| indicates 
summation over j = 1, 2, 3. 

Torsion is a rank three tensor antisymmetric in the last two indices. Because 
we are working in dimension three, it is convenient, as in [3], to apply the Hodge 
star in the last two indices and deal with the rank two tensor 

T a := \ T^ s e lS0 ^det^ (3.14) 

instead. Substituting p. 121) into (|3.14l) we get 

T = ej® curie 5 ', (3.15) 

where 

(cm\e 3 )p := (^de 3 )^ = - (de 3 ) 15 e 1&p y/detg^ . (3.16) 

4 Relation between curvature of the U(l) connec- 
tion and torsion of the teleparallel connection 

This section is devoted to the examination of the integrand in formula (|2.7|) . 
Recall that the curly brackets in this integrand denote the Poisson bracket on 
matrix-functions (|2.5jl . 

As explained in Section 5 of [10], the expression — is the scalar 
curvature of the U(l) connection generated by the eigenspace spanu + of the 
principal symbol. This curvature term appears in the general setting of a first 
order elliptic system. A feature of the particular case (jl.9|) - (11.12p considered 
in the current paper is that the scalar curvature of the U(l) connection can 
be expressed via torsion of the teleparallel connection. This is a substantial 
simplification. The teleparallel connection is a simpler geometric object than 
the U(l) connection because the coefficients of the teleparallel connection do 
not depend on the dual variable (momentum), i.e. they are "functions" on the 
base manifold M. The relationship between the two connections is established 
by the following lemma. 

Lemma 4.1 The scalar curvature of the U(l) connection is expressed via the 
torsion of the teleparallel connection, metric and topological charge as 

-WlhO-j ^ . («) 
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Recall that the topological charge c = ±1 is defined in accordance with 
formula (|1.17p or, equivalently, in accordance with formula (|3.5p . 

Proof of Lemma \4-l\ We give the proof for the case 

c = +l. (4.2) 

There is no need to give a separate proof for the case c = — 1 as the two cases 
reduce to one another by means of a) the observation that torsion (|3.12p is 
invariant under inversion of the frame and b) the identity 

{[v + ]*,v+} + {[v-]*,v-} = 0, (4.3) 

where v~ [x, £) is the normalized eigenvector of the principal symbol correspond- 
ing to the negative eigenvalue. Formula (|4.3p is a special case of formula (1.22) 
from [10] . 

We fix an arbitrary point Q € T'M and prove formula (|4.ip at this point. 
As the LHS and RHS of (|4.1I) are invariant under changes of local coordinates x, 
it is sufficient to prove formula (|4.ip in Riemann normal coordinates, i.e. local 
coordinates such that x — corresponds to the projection of the point Q onto 
the base manifold, g Ml /(0) = 8^ and ^pr(O) = 0. Moreover, as the formula we 
are proving involves only first partial derivatives in x, we may assume, without 
loss of generality, that 

g^{x) = <V (4.4) 

for all x in some neighbourhood of the origin. In other words, it is sufficient to 
prove formula (|4.ip for the case of Euclidean metric. 

As both the LHS and RHS of (|4.ip have the same degree of homogeneity 
in £, namely, —1, it is sufficient to prove formula (|4.ip for £ of norm 1. Moreover, 
by rotating our Cartesian coordinate system we can reduce the case of general 
£ of norm 1 to the case 

^=(0 1). (4.5) 

There is one further simplification that can be made: we claim that it is 
sufficient to prove formula (|4.ip for the case when 

e j a (0) = 6 j a , (4.6) 

i.e. for the case when at the point x — the elements of the frame are aligned 
with the coordinate axes. This claim follows from the observation that the 
LHS of formula (|4.ip is invariant under rigid special unitary transformations of 
the column- function v + . v + i— > Rv + , where "rigid" refers to the fact that the 
matrix R G SU(2) is constant. Of course, the column-function Rv + is no longer 
an eigenvector of the original principal symbol, but a new principal symbol 
obtained from the old one by the rigid special orthogonal transformation of the 
frame (|A.14[) with the 3x3 special orthogonal matrix O expressed in terms 
of the 2x2 special unitary matrix R in accordance with (|A.15p . One can 
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always choose the special unitary matrix R so that at the point x — the 
elements of the new frame are aligned with the coordinate axes (in fact, there 
are two possible choices of R which differ by sign) . It remains only to note that 

direct inspection of formula (I3.12[) shows that torsion is also invariant under 

* 

rigid special orthogonal transformations of the frame, and, hence, the tensor T 
defined by formula (|3.14[) and appearing in the RHS of formula (|4. 1 [) is invariant 
under rigid special orthogonal transformations of the frame as well. 

Having made the simplifying assumptions (|4.4I) - (|4.6|) . we are now in a posi- 
tion to prove formula (|4.ip . 

Let us calculate the RHS of (|4.1[) first. In view of (|4.6[) we have, in the linear 
approximation in x, 

'e^ix) e x 2 (x) ei 3 (x)\ / 1 w 3 (x) -w 2 {x)\ 

e2 1 (x) e 2 2 (x) e 2 3 (x) = -w 3 (x) 1 w^x) , (4.7) 

K e 3 1 (x) e 3 2 (x) e 3 3 (x)J \w 2 (x) -w 1 (x) 1 / 

where w is some smooth vector- function which vanishes at x = 0. Formula (|4.7p 
is the standard formula for the linearization of an orthogonal matrix about the 
identity; see also formula (10.1) in [3]. Note that in Cosserat elasticity literature 
the vector-function w is called the vector of microrotations. Substituting (|4.7|) 
into (|3.15p and (|3.16p we get, at x = 0, 

T aP = dwp /dx a - 5 af} div w, (4.8) 

which is formula (10.5) from [3]. Here we freely lower and raise tensor indices 
using the fact that the metric is Euclidean (in the Euclidean case (|4.4|) it does 
not matter whether a tensor index comes as a subscript or a superscript). Sub- 
stituting (|3~8]) and ([4~5]) into the RHS of (|4~H we get, at our point Q E T'M, 

li^^^-¥ dwl/dxl+dw2/dx2) - (4 - 9) 

Let us now calculate the LHS of (|4.ip . The equation for the eigenvector 
v + (x,{;) of the principal symbol is 



e 3 Q £aH!£ll (ei-ie 2 ) Q £ Q \ fv+ 



(4.10) 



In view of (|4.5p and (|4.6p the (normalized) solution of (|4.10p at our point 
0, 

v + ^e % f'v + , (4.11) 



Q G T'M is v + • •■ ( n ) • 01' course, our r + (,r.i) is < It •fined up t.» I he {•an.nc 
transformation 



where 

(j)+:T'M^R (4.12) 



13 



is an arbitrary smooth function, however the LHS of (|4.1[) is invariant under 
this gauge transformation. We now perturb equation (|4.10[) about the point 
Q G T'M, that is, about x = 0, £ = (0 l), making use of formula ([iTfjl . 
which gives us the following equation for the increment 5v + of the eigenvector 
v + (x,{;) of the principal symbol: 

0WS«A / -w 2 (x) -iw 1 (x)\ fl 

-2 J \Sv+) + \-w 2 (x) +iw 1 (x) ) [O 

6£i - i6&\ A 



<*6 + -256 / V 



0. 



or, equivalently, 



<5v+ = ^(-w 2 (x) + iw 1 {x) + 56 + itffa). (4.13) 



Formula (|4. 13[) has to be supplemented by the normalization condition 
|i> + (a;,6l = 1, which in its linearized form reads 

Re5v+=0. (4.14) 

Formulae (|4.14[) and (|4.13[) define Sv + modulo an arbitrary ImSv^, with this 
degree of freedom being associated with the gauge transformation (|4. 1 1|) , (|4. 12|) . 
Without loss of generality we may assume that the gauge is chosen so that 

Imck>+=0. (4.15) 

Combining formulae (|4.14[) . (14.151) and (|4.13p we get 

5V+ = 2 {-w 2 (x) + iw\x) + (56 + W&) ' (4 ' 16) 

Recall that the w appearing in this formula is some smooth vector-function 
which vanishes at x = 0. 

Differentiation of (|4.16p gives us 



dv+ _ 1 / 

dx°< ~ 2 \-dw 2 /dx a + idw 1 /dx c 

dv+ 1 f0\ dv+ 1 f0\ dv 



96 2 \ij ' 96 2 Vv ' 96 

Formulae (|4TT7|) and P~T5]) imply that at our point Q e T'M 



(4.17) 
(4.18) 



-i{[v+]*,v+} = -^{dw 1 /dx 1 +dw 2 /dx 2 ). (4.19) 

Comparing formulae (|4.9[) and (|4.19[) and recalling (|4.2p . we arrive at the 
required result (|4.ip . □ 
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5 Integration of the curvature term 

Substituting into (|2T7| we get 

h+(a,0<l 

Recall that ft, + (a;,£) is given by formula (|1.14p . 

The tensor T can be decomposed into pure trace and trace-free pieces, i.e. 

T a/3 = -g^h^ + ( T ap - -g af3 h 1 J . (5.2) 



It is easy to see that the trace- free piece (second term in the RHS of (|5.2p ) does 
not contribute to the integral in (|5.ip . hence formula (|5.1[) becomes 

62(x) = ^TT 7 (aj) J d£ = ^{b iy /M^)(x). (5.3) 

But formulae fOD]) . (|3T4j) and (f3~T4|) imply that 

T\ = 3 * T ax . (5.4) 

Combining formulae (|2T2|) . ([23]) . (|S3| and (|5T4|) we arrive at formula (fl~19|) . 
This completes the proof of Theorem 11.11 



6 The subprincipal symbol of the massless Dirac 
operator 

In this section we calculate the subprincipal symbol of the massless Dirac opera- 
tor, which prepares the ground for the proof of Theorem 1 1.2 1 in the next section. 
In view of Remark 2.1.10 from [22J, defining the subprincipal symbol for the 
massless Dirac operator on spinors (|A.3j) is problematic, hence, we work with 
the massless Dirac operator on half-densities (|A.19j) . For the sake of brevity we 
denote the massless Dirac operator on half-densities by A rather than by Wi/ 2 ■ 

Lemma 6.1 The subprincipal symbol of the massless Dirac operator on half- 
densities tA.19\) is 

A auh (x) = ^(*T°*(x))l, (6.1) 

where c = ±1 is the topological charge VJ.5\) , T ax is axial torsion il.20\) , * is 
the Hodge star \3.$ and I is the 2x2 identity matrix. 



15 



Proof We give the proof of (|6.ip for the case (|4.2[) . There is no need to give 
a separate proof for the case c = — 1 as the two cases reduce to one another by 
inversion of the frame: the full symbol of the massless Dirac operator on half- 
densities changes sign under inversion of the frame and hence its subprincipal 
symbol changes sign under inversion of the frame, whereas torsion ()3.12p is 
invariant under inversion of the frame. 

According to formula (1.2) from [12J the subprincipal symbol is defined as 

A suh := A + ^(Ai) x «. €o , (6.2) 

where A\(x,£) and Ao(x) are the homogeneous (in £) components of the full 
symbol A(x,£) — Ai(x,£) + Aq(x) of our first order differential operator, with 
the subscript indicating degree of homogeneity. For the massless Dirac operator 
on half-densities (|A.19|) these homogeneous components read (|3.ip and 

respectively. Note that in writing down (|6.3p we used the standard formula 

1 d det g^ v _ { P \ 
2det 5ftA dx a ~\ap]' 

Our task is to substitute (|3. 1|) and (|6.3p into (|6.2p . 

We fix an arbitrary point P 6 M and prove formula (|6.ip at this point. As 
the LHS and RHS of (|6.ip are invariant under changes of local coordinates x, it 
is sufficient to check the identity (|6.ip in Riemann normal coordinates, i.e. local 
coordinates such that x = corresponds to the point P, g^ v (0) = 8^ v and 
|^pr(0) = 0. Moreover, as the principal symbol is linear in £ and the formula we 
are proving involves only first partial derivatives in x, we may assume, without 
loss of generality, that we have (|4.4p for all x in some neighbourhood of the 
origin. In other words, it is sufficient to prove formula (|6.ip for the case of 
Euclidean metric. Furthermore, by rotating our Cartesian coordinate system we 
can achieve (|4.6I) . which opens the way to the use, in the linear approximation 
in x, of formula (|4.7p . 

Substituting (|4.7p into (jA.ip . we get, in the linear approximation in x, 

1 ( w 2 1 + iw 3 \ 
° = [l iw 3 -u? ) = ai ' 

2 / —w 1 —i + w 3 
\i + w w 

q / 1 — iw 1 — w 2 \ ,„ .s 

v = ■ 1 2 1 = 0-3 . 6.4 

\iw —w -1 J 

Recall that the w appearing in this formula is some smooth vector-function 
which vanishes at x = 0. 
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Substitution of (|6.4[) into (|3.ip and (|fc>. 3[) gives us 



Ai(a;,0 



6 6 - ifa 

w 2 £i - w; 1 ^ iw 3 £i + w 3 £ 2 + (-iw 1 - w 2 )£ 3 

iw 3 ^ + u> 3 £ 2 + (iw 1 - u> 2 )£ 3 -w 2 £i + w 1 ^ 



(6.5) 



i fO i\ fO l\ / dw 2 /dx 1 idw 3 /dx 1 \ . . 

A °^- l\i o)[i 0j{-idw 3 /dx 1 -dw 2 /dx 1 ) + '--- {b - b > 

Here formula (|6.5p is written in the linear approximation in x, whereas formula 
(|6.6p displays, for the sake of brevity, only one term out of nine (the one cor- 
responding to a = j3 = 1 in (|6.3p ). with the remaining eight terms concealed 
within the dots .... Note also that the Christoffel symbols disappeared because 
of our assumption that the metric is Euclidean. 
Substituting ([676]) and (|675|) into lj6T2]t . we get 



A sub (0) = -i(divw)7. (6.7) 

But, according to (|4.8|) . 

T 7 7 (0) = -2divu;. (6.8) 
Formulae ([577]) . (JSTSJ, |57J| and (O imply formula {BHJ at x = 0. □ 



7 Proof of Theorem 11.21 

As Theorem II .21 is an if and only if theorem, our proof comes in two parts. 



Part 1 of the proof Let A be a massless Dirac operator on half-densities. 
We need to prove that a) the subprincipal symbol of this operator, A su \ ) (x), is 
proportional to the identity matrix and b) the second asymptotic coefficient of 
the spectral function, b(x), is zero. The required result follows from Lemma |6. II 
and Theorem ll.il 



Part 2 of the proof Let A be a differential operator such that a) the sub- 
principal symbol of this operator, A su ^(x), is proportional to the identity matrix 
and b) the second asymptotic coefficient of the spectral function, b(x), is zero. 
We need to prove that A is a massless Dirac operator on half-densities. 

Theorem 1 1.1 1 implies that the subprincipal symbol of our operator A is given 
by formula (|6.ip . Let ej be the frame corresponding to the principal symbol of 
the operator A, see formulae (|3.ip and (|3.2p . Now, let B be the massless Dirac 
operator on half-densities corresponding to the same frame. Then the principal 
symbols of the operators A and B coincide. But Lemma 16.11 implies that the 
subprincipal symbols of the operators A and B coincide as well. A first order 
differential operator is determined by its principal and subprincipal symbols, 
hence, A = B. □ 
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8 Explicit formula for axial torsion 



Torsion is a rank three tensor antisymmetric in the last two indices. It is known 
[51 ITS] that torsion has three irreducible pieces. Only one of the three irreducible 
pieces of torsion, namely, the piece which theoretical physicists label by the 
adjective "axial", appears in our spectral theoretic results, see Theorem 11.11 and 
Lemma |6. 11 It is also interesting that axial torsion is the irreducible piece which 
is used when one models the massless neutrino [11, or the electron [B] by means 
of Cosserat elasticity. 

Axial torsion is defined as the totally antisymmetric piece of the torsion 
tensor, see formula (|1.20|) . This means that axial torsion is a 3-form. In view of 
the importance of axial torsion, we give an explicit formula for its Hodge dual in 
terms of the principal symbol A%(x, £). Formulae (|3.15p . (|3.16p and (15. 4p imply 



o 

- e\ de l 2 /dx 3 - e k 2 de^/dx 1 - e k 3 de\/dx 2 ] . (8.1) 

Here the coframe e k is determined from the principal symbol in accordance with 
formulae (|3.ip . (|3.2p and (|3.9p , whereas the contravariant metric tensor <? Q/3 is 
determined from the principal symbol in accordance with formula (| 1 . 1 3[) . 
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A The massless Dirac operator 

Let M be a 3-dimensional connected compact oriented manifold equipped with 
a Riemannian metric g a p, a, (3 = 1,2,3 being the tensor indices. Note that we 
are more prescriptive in this appendix than in the main text of the paper: in the 
main text orientability emerged as a consequence of the existence of a principal 
symbol and the metric was defined via the principal symbol, whereas in this 
appendix orientability and metric are introduced a priori. 

We work only in local coordinates with prescribed orientation. 

It is known [231 [T5] that a 3-dimensional oriented manifold is paralleliz- 
able, i.e. there exist smooth real vector fields ej, j = 1,2,3, that are linearly 
independent at every point x of the manifold. (This fact is often referred to 
as Steenrod's theorem.) Each vector ej(x) has coordinate components ej a (x), 
a = 1,2,3. Note that we use the Latin letter j for enumerating the vector 
fields (this is an anholonomic or frame index) and the Greek letter a for enu- 
merating their components (this is a holonomic or tensor index). The triple 
of linearly independent vector fields ej, j — 1, 2, 3, is called a frame. Without 
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loss of generality we assume further on that the vector fields ej are orthonor- 
mal with respect to our metric: this can always be achieved by means of the 
Gram-Schmidt process. 



Define Pauli matrices 



(j a {x) :=s J e 3 a (x), (A.l) 



where 



•'-(! ;)-- M° 7)-' Mi -i)="- <-> 

In formula (|A.1|) summation is carried out over the repeated frame index j — 
1,2,3, and a — 1,2,3 is the free tensor index. 

The massless Dirac operator is the matrix operator 

— (^+H^ + {»)' <"» 

where summation is carried out over a, /?, 7 = 1,2,3, and 



Z 3 1 . = i ^5 /" ^SV , d 9as _ dg a 



07 



2 V dx a dx~< dx s 



(A.4) 



are the Christoffel symbols. Here and throughout this appendix we raise and 
lower tensor indices using the metric. Note that we chose the letter "W" for 
denoting the massless Dirac operator because in theoretical physics literature it 
is often referred to as the Weyl operator. 

Formula (|A.3[) is the formula from |11| . only written in matrix notation 
(i.e. without spinor indices). Note that in the process of transcribing formulae 
from |llj into matrix notation we used the identity 

ea a e = (a a ) T , (A.5) 

a = 1, 2, 3, where 

; - 1 ) <A,) 

is the "metric spinor". The identity (|A.5[) gives a simple way of raising/lowering 
spinor indices in Pauli matrices in the non-relativistic (a ^ 0) setting. 

Our definition (IA.3P of the massless Dirac operator is a special case of the 
definition from |13| . The two definitions coincide when we work with a Spin 
connection as opposed to a Spin c connection, see Propositions 2.14 and 2.15 in 
[T3] for details. 

Throughout this paper we work in dimension 3. The definition of the mass- 
less Dirac operator acting over a Riemannian manifold of arbitrary dimension 
can be found, for example, in [THl HH ITS] . 

Physically, our massless Dirac operator (|A.3[) describes a single massless 
neutrino living in a 3-dimensional compact universe M. The eigenvalues of the 
massless Dirac operator are the energy levels. 
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The massless Dirac operator (|A.3[) acts on 2-columns v — (v\ v%\ of 
complex- valued scalar functions. In differential geometry this object is referred 
to as a (Weyl) spinor so as to emphasize the fact that v transforms in a particular 
way under transformations of the orthonormal frame ej. However, as in our 
exposition the frame ej is assumed to be chosen a priori, we can treat the 
components of the spinor as scalars. This issue will be revisited below when we 
state Property 4 of the massless Dirac operator. 

We now list the main properties of the massless Dirac operator. We state 
these without proofs. The proofs can be found in Appendix 3. A of [5] or in |13) . 

Property 1. The massless Dirac operator is invariant under changes of 
local coordinates x, i.e. it maps 2-columns of smooth scalar functions M C 2 
to 2-columns of smooth scalar functions M — > C 2 regardless of the choice of 
local coordinates. 



Property 2. The massless Dirac operator is formally self-adjoint (symmet- 
ric) with respect to the inner product 

w*v det g a /3 dx (A. 7) 

on 2-columns of smooth scalar functions v,w : M — > C 2 . 

Property 3. The massless Dirac operator W commutes 

C(Wv) = WC(v) (A.8) 

with the antilinear map 

v ^ C(v) := ev, (A.9) 

where e is the "metric spinor" (|A.6I) . In theoretical physics the transformation 
(|A.9[) is referred to as charge conjugation [H I13| . 

Formula (|A.8|) implies that v is an eigenfunction of the massless Dirac op- 
erator corresponding to an eigenvalue A if and only if C(i>) is an eigenfunction 
of the massless Dirac operator corresponding to the same eigenvalue A. Hence, 
all eigenvalues of the massless Dirac operator have even multiplicity. Moreover, 
any eigenfunction v and its "partner" C(v) make the same contribution to the 
spectral function (|1.3[) at every point x of the manifold M. 

If, as in |13j . we introduce a magnetic field, then we lose the commutation 
property (|A.8|) and the double eigenvalues split up. This indicates that the 
double eigenvalues of the massless Dirac operator correspond to the two different 
spins. 

Property 4. This property has to do with a particular behaviour under 
SU(2) transformations. Let R : M — >• SU(2) be an arbitrary smooth special 
unitary matrix- function. Let us introduce new Pauli matrices 

a a := Ra a R* (A.10) 
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and a new operator W obtained by replacing the a in (|A.3I) by a. It turns out 
(and this is Property 4) that the two operators, W and W, are related in exactly 
the same way as the Pauli matrices, a and a, that is, 

W = RWR*. (A.ll) 

We now examine the geometric meaning of the transformation (|A.10[) . Let 
us expand the new Pauli matrices a with respect to the basis (|A.2[) : 

d a {x) = s j ij a {x). (A.12) 

Formulae (|A.1[) , (|A.12[) and (|A.10[) give us the following identity relating the 
new vector fields ej and the old vector fields ej : 

Rs k R*e k = s j e, . (A.13) 

Resolving (|A.13[) for ej we get 

< s ()/■'!.-. (A.14) 
where the real scalars Ojk are given by the formula 

0/ = ±tT( Sj Rs k R*). (A.15) 

Note that in writing formulae (|A.13j) and (|A.14j) we chose to hide the tensor 
index, i.e. we chose to hide the coordinate components of our vector fields. Say, 
formula (|A.14[) written in more detailed form reads e.j a = Oj k ek a ■ 

The scalars (|A.15[) can be viewed as elements of a real 3x3 matrix-function 
O with the first index, j, enumerating rows and the second, k, enumerating 
columns. It is easy to check that this matrix-function O is special orthogonal. 
Hence, the new vector fields ij are orthonormal and have the same orientation as 
the old vector fields ej . We have shown that the transformation (|A.10[) has the 
geometric meaning of switching from our original oriented orthonormal frame 
ej to a new oriented orthonormal frame ej . 

Formula (|A.15I) means that the special unitary matrix R is, effectively, a 
square root of the special orthogonal matrix O. It is easy to see that for a given 
matrix O £ S0(3) formula (|A.15[) defines the matrix R 6 SU(2) uniquely up 
to sign. This observation allows us to view the issue of the geometric meaning 
of the transformation (|A.10[) the other way round: given a pair of orthonormal 
frames, ej and ej, with the same orientation (i.e. with sgndet ej a — sgndet ej a ), 
we can recover the special orthogonal matrix- function 0(x) from formula (|A.14p 
and then attempt finding a smooth special unitary matrix-function R(x) satis- 
fying (|A.15[) . Unfortunately, this may not always be possible due to topological 
obstructions. We can only guarantee the absence of topological obstructions 
when the two frames, ej and ej, are sufficiently close to each other, which is 
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equivalent to saying that we can only guarantee the absence of topological ob- 
structions when the special orthogonal matrix-function 0(x) is sufficiently close 
to the identity matrix for all x £ M . 

We illustrate the possibility of a topological obstruction by means of an 
explicit example. Consider the unit torus T 3 parameterized by cyclic coordinates 
x a , a = 1,2,3, of period 2ir. The metric is assumed to be Euclidean. Define a 
pair of orthonormal frames 

and 



:= S 3 a (A.16) 



^cosfc3ir 3 \ /— sinA^x 3 ^ 

e x a := I sink 3 x 3 , e 2 a := cosk^x 3 I , e 3 a := I I , (A.17) 



where £3 is an odd integer. Let W and W be the massless Dirac operators 
corresponding to the frames (|A.16I) and (|A.17|) respectively. We claim that 
there does not exist a smooth matrix-function R : T 3 — > SU(2) which would 
give (|A.15[) . where 0(x) is the special orthogonal matrix- function defined by 
formula (|A.14j) . We justify this claim in two different ways. 

Justification 1. Resolving the system (|A.14j) - (|A.17jl locally for R, we get 

^ 3 H±( e T 3 jj)> (A.18) 

and this solution is unique modulo choice of sign; here the freedom in the choice 
of sign is not surprising as SU(2) is the double cover of SO (3). Formula (|A.18[) 
defines a continuous single- valued matrix-function on the unit torus T 3 if and 
only if the integer k 3 is even, which it is not. 

Justification 2. It is sufficient to show that the two operators, W and W, 
have different spectra. Straightforward separation of variables shows that zero 
is an eigenvalue of the operator W but not an eigenvalue of the operator W. 

One can generalize the above example by introducing rotations in three 
different directions, which leads to eight genuinely distinct parallelizations. See 
also [24] page 524 or [2] page 21. 

Let us emphasize that the topological obstructions we were discussing have 
nothing to do with Stiefel-Whitney classes. We are working on a parallelizable 
manifold and the Stiefel-Whitney class of such a manifold is trivial. The topo- 
logical issue at hand is that our parallelizable manifold may be equipped with 
different spin structures. 

We say that two massless Dirac operators, W and W, are equivalent if there 
exists a smooth matrix-function R : M — > SU(2) such that the corresponding 
Pauli matrices, a a and a a , are related in accordance with (|A.10[) . In view of 
Property 4 (see formula ([A. lip ) all massless Dirac operators from the same 
equivalence class generate the same spectral function (jl.3[) and the same count- 
ing function (|1.4[) . so for the purposes of our paper viewing such operators as 
equivalent is most natural. 
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As explained above, there may be many distinct equivalence classes of mass- 
less Dirac operators, the difference between which is topological. Studying the 
spectral theoretic implications of these topological differences is beyond the 
scope of our paper. The two-term asymptotics (ll.22[) and (|1.23p derived in the 
main text of our paper do not feel this topology. 

In theoretical physics the SU(2) freedom involved in defining the massless 
Dirac operator is interpreted as a gauge degree of freedom. We do not adopt 
this point of view (at least explicitly) in order to fit the massless Dirac operator 
into the standard spectral theoretic framework. 

We defined the massless Dirac operator (|A.3[) as an operator acting on 2- 
columns of scalar functions, i.e. on 2-columns of quantities which do not change 
under changes of local coordinates. This necessitated the introduction of the 
density */ det g a p in the formula (|A.7[) for the inner product. In spectral theory 
it is more common to work with half-densities. Hence, we introduce the operator 



which maps half-densities to half-densities. We call the operator (|A.19[) the 

massless Dirac operator on half- densities. 

B The spectrum for the torus and the sphere 

In this appendix we examine the massless Dirac operator on the unit torus 
T 3 and the unit sphere § 3 and compare our asymptotic formulae (|1.22[) and 
(jl.23p with known explicit formulae. The torus is assumed to be equipped with 
Euclidean metric (see also Appendix whereas the sphere is assumed to be 
equipped with metric induced by the natural embedding of § 3 in Euclidean space 
M 4 . Note that in view of the obvious symmetries of the torus and the sphere 
the scalar function e(A, x, x)j det g a p (x) is constant (see also Remark ll.2|) . so 
formulae (|1.22j) and (|1.23j) are in this case equivalent, in the sense that they 
follow from one another. Hence, we will be dealing with formula Q1.23P only. 
We have VolT 3 = (27r) 3 , so for the torus formula (|1.23p reads 



The nonperiodicity condition (see Definitions 8.3 and 8.4 in [10J ) is fulfilled for 
the torus, so, according to Theorem 8.4 from [TU], the asymptotic formula (|B.1[) 
holds as it is, without mollification. 

In order to test formula (jB.ip we calculate the spectrum of the massless Dirac 
operator on T 3 explicitly. We do this first for the spin structure associated with 
the frame (|A.16|) . Then the spectrum is as follows. 

• Zero is an eigenvalue of multiplicity two. 

• For each m E 1? \ {0} we have the eigenvalue ||m|| and unique (up to 
rescaling) eigenfunction, with eigenfunctions corresponding to different m 
being linearly independent. 



W 1/2 := (det g^Y^W (det g^y 1 / 4 



(A.19) 




ttA 3 + o(A 2 ). 



(B.l) 
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• For each m € Z 3 \ {0} we have the eigenvalue — |jm|| and unique (up to 
rescaling) eigenfunction, with eigenfunctions corresponding to different m 
being linearly independent. 

Hence, ^V(A) + 1 is the number of integer lattice points inside a 2-sphere of 
radius A in R 3 centred at the origin. According to [T7] the latter admits the 
asymptotic expansion 

^A 3 + O e (A 21 / 16+E ) (B.2) 

as A — > +00, with e being an arbitrary positive number. This agrees with our 
asymptotic formula (IB.1[) . 

Let us now consider the spin structure associated with the frame (| A. 17|) . 
Then the spectrum is as follows. 

• For each m £ Z 3 we have the eigenvalue ||m — (0,0, 1/2) || and unique (up 
to rescaling) eigenfunction, with eigenfunctions corresponding to different 
m being linearly independent. 

• For each m S Z 3 we have the eigenvalue — \\m— (0, 0, 1 /2) || and unique (up 
to rescaling) eigenfunction, with eigenfunctions corresponding to different 
m being linearly independent. 

Hence, N(X) is the number of integer lattice points inside a 2-sphere of radius 
A in M 3 centred at (0,0,1/2). Here the sphere is shifted from the origin so 
one cannot apply the result from |17] . However, as the shift is rational, one 
can reduce the problem to counting integer lattice points in a rational ellipsoid 
centred at the origin, and an application of the result from [8] gives us for the 
shifted sphere the same asymptotic expansion (|B.2j) as for the sphere centred at 
the origin. 

As explained in Appendix [A] the unit torus T 3 admits a total of eight dif- 
ferent spin structures. For each of these the problem of counting positive eigen- 
values of the massless Dirac operator reduces to counting integer lattice points 
inside a 2-sphere of radius A in M 3 (possibly, shifted from the origin by a ra- 
tional shift), so in all eight cases we do get (|B.1|) . In fact, we can replace the 
remainder o(A 2 ) in (IB.1[) by O e (A 21 / 16+e ) and this holds for all eight different 
spin structures. 

In the remainder of this appendix we examine the massless Dirac operator 
on the unit sphere § 3 . We have Vol§ 3 = 27r 2 , so for the sphere formula (|1.23|) 
reads 

^(A) = y+o(A 2 ). (B.3) 

The nonperiodicity condition fails for the sphere because all geodesies are closed 
with period 2ir, so formula (|B.3[) cannot be used in its original form and has to 
be mollified, see Remark 1 1.1 1 We will deal with the mollification issue later and 
give explicit formulae for the eigenvalues first. 

It is known that § 3 admits a unique spin structure, see Section 5 in [2]. 
The spectrum of the massless Dirac operator on S 3 has been computed by 
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different authors using different methods [26 1 [27] HI H] ancl reads as follows: the 
eigenvalues are 

±(k+~y fc = l,2,..., (B.4) 

and their multiplicity is 

k(k + l). (B.5) 

The mollification procedure from Section 7 of [ID] goes as follows. Put 
N(X) :— for A < and take an arbitrary real-valued even function p(A) from 
Schwartz space 5(R) whose Fourier transform p(t) satisfies conditions p(0) = 1 
and supp p C (—2tt, 2tt). Then, according to Theorem 7.2 from [TDJ, the mollified 
version of formula (|B.3|) reads 

J N(\- f i)p(p)dfi=^- + 0(X) 

and this result holds notwithstanding the failure of the nonperiodicity condition. 
However, for the sphere there is a much simpler way of testing our asymptotic 
formula. Let A > 2 be integer. Taking an integer A puts us exactly in the middle 
of the gap between two consecutive clusters of eigenvalues, see formulae (|B.4|) 
and (|B.5[) , and achieves the same averaging effect as convolution with a function 
from Schwartz space. For integer A > 2 we get 

N(\) = '£k(k + 1) = --- 

k=l 

which agrees with our asymptotic formula (IB.3[) . 
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